co = Crystal.create_hematite_rhombohedral;
co.convert_rhombohedral_to_hexagonal;
%co.replicate([3 3 3]);

clf;
co.plot_atoms;

co.write_xyz_file('hematite');


pos = POSCAR('hematite', 1, co.LatticeVectors,...
        co.NumAtoms, 0, '',...
        co.Positions)

pos.write;

%co2 = Crystal.create_hematite_rhombohedral;
%pos2 = POSCAR('hematite', 1, co2.LatticeVectors,...
%        co2.NumAtoms, 0, '',...
%        positionCoordinates, initialVelocities, velocityCoordinates)
